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Interacting many-particle systems with a mean-field one body part plus a chaos generating random 
two-body interaction having strength A, exhibit Poisson to GOE and Breit-Wigner (BW) to Gaussian 
transitions in level fluctuations and strength functions with transition points marked by A = A c 
and \ = Xp, respectively; Xp » A c . For these systems theory for matrix elements of one-body 
■ transition operators is available, as valid in the Gaussian domain, with A > Af, in terms of orbitals 

' occupation numbers, level densities and an integral involving a bivariate Gaussian in the initial and 

. final energies. Here we show that, using bivariate ^distribution, the theory extends below from the 

Gaussian regime to the BW regime up to A = A c . This is well tested in numerical calculations for 
six spinless fermions in twelve single particle states. 
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Two-body random matrix ensembles apply in a generic way to finite interacting many fermion systems such as nuclei 
[1, 2], atoms [3, 4], quantum dots [5], small metallic grains [6] etc. A common feature of all these systems is that their 
^""^ ' hamiltonian (H) consists of a mean-field one-body [h(l)J plus a complexity generating two-body [V(2)] interaction. 
With this, one has EGOE(l+2), the embedded Gaussian orthogonal ensemble of one plus two-body interactions 
operating in many particle spaces [2]; for more complete definition of EGOE(l+2) for m fermions in N single particle 
states see [7]. Most significant aspect of EGOE(l+2) is that as A, the strength of the random (represented by GOE) 
i—i two-body interaction [in H = h(l) + XV (2)], changes, in terms of state density, level fluctuations, strength functions 
and entropy [8], the ensemble admits three chaos markers. Firstly, it is well known that the state densities take 
Gaussian form, for large enough to, for all A values [9]. With A increasing, there is a chaos marker A c such that for 
y—i " A > A c the level fluctuations follow GOE, i.e. A c marks the transition in the nearest neighbor spacing distribution 
from Poisson to Wigner form [10]. As A increases further from A c , the strength functions (for h(l) basis states) change 
from Breit-Wigner (BW) to Gaussian form and the transition point is denoted by Af [11]. The A c < A < Af region 
is called BW domain and the A > Af region is called Gaussian domain. As we increase A much beyond Af, there 
is a chaos marker At around which different definitions of entropy, temperature etc. will coincide and also strength 
functions in h(l) and V(2) basis will coincide. Thus A ~ At region is called the thermodynamic region [12, 13]. 
• With the three chaos markers A c , Af and At, EGOE generates statistical spectroscopy, i.e. smoothed forms for 
state densities, orbit occupancies, strength sums [for example Gamow- Teller (GT) sums in nuclei, electric dipole (El) 
sums in atoms], transition strengths themselves [for example: electric quadrupole(i?2), magnetic dipole (Ml) and GT 
strengths in nuclei, El strengths in atoms and molecules etc.], information entropy in wavefunctions and transition 
strengths etc. The EGOE Gaussian state densities are being used to generate a theory (valid for A > A c ) for level 
densities with interactions [14]. Similarly, theory for orbit occupancies and strength sums, as valid in BW to Gaussian 
regimes (i.e. for A > A c ) has been developed [13]. However, for transition strengths (experimentally they are most 
important for probing wavefunctions structure of a quantum system) , a theory valid only in the Gaussian domain is 
available [15-17]. Although a theory was given by Flambaum et al for BW domain [3, 18, 19], it is well known to 
underestimate the exact values by a factor of 2 [17, 18]. Thus, a major gap (see the discussion in [17]) in understanding 
transition strengths is in extending the theory that works in the Gaussian domain, well into the BW domain. The 
purpose of this paper is to show that the bivariate i-distribution known in statistics will bridge this gap. As in Refs. 
[17, 19], we restrict ourselves to one-body transition operators. 

Given a Hamiltonian H and its m-particle eigenstates \E), the transition strengths generated by a one-body 
transition operator O are denoted by \{Ef \ O | Ei)\ ; O = J2 a ~ eapa^ap where a^ a creates a particle in the single 
particle state a and ap destroys a particle in the state (3. Now the bivariate strength densities are defined by 

I&7o( E *> E f) = i(°' S ( H E f)°^ H E *))T = ((0"0)) m P^h(E t ,E f ) . (1) 

In Eq. (1), (( )) denotes trace. Note that Iuv,o is square of the matrix elements of O in H eigenstates 
weighted by the state densities at the initial and final energies and the corresponding puv\0 is normalized to unity. 
Moreover, one-body transition operators O will not change m. The bivariate moments of puv,o are defined by 
M m = {(O^HiOHP))" 1 1 {{O^O))" 1 . With Mm = e t and M m = e f defining the centroids of its two marginals, the 
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bivariate central moments of pu V \0 are given by 

p pq = «0t (H - e f )« 0{H- e t r)) m / ((0^0)) m . (2) 

Most important of these are = of and ^02 = O/, the variances of the two marginals and ( = pn/criif, the 
bivariate correlation coefficient. 

For EGOE(l+2), going well into the Gaussian domain [then EGOE(l+2) will be effectively EGOE(2)], it is well 
established that the bivariate strength densities take bivariate Gaussian form (this applies to nuclei [15, 16]), 



Pbiv,o(Ei, Ef) pbi v -g ; o(Ei,Ef,ei,ef,ai,af,Q 




(3) 
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An immediate question is how to extend this result well into the BW domain and up to A c (note that GOE fluctuations 
operate for A > A c and hence in this regime it is possible to consider smoothed transition strengths). In a recent 
work, Angom et al [20] showed that strength functions covering the BW to Gaussian regimes can be well represented 
by Student's t-distribution. Following this result, here we conjecture that the bivariate strength density puv.p in Eq. 
(1) can be represented by the bivariate i-distribution, 

pbiv-t;o(Ei,Ef;e i ,ef,ai,(T2,C, l/ ) = - 1 t. , 2 x 

i+^{(^) 2 -^)(^n^r rii? (4) 



V > 1 . 



Properties of puv-t are given in [21, 22]. Most important is that for (y = 1, £ = 0), puv-t gives bivariate BW (called 
bivariate Cauchy in statistics) distribution and as v — > 00, puv-t becomes bivariate Gaussian. Thus it has the correct 
limiting forms and the intermediate shapes are largely determined by the v parameter. The marginal distributions 
of puv-t are easily seen to be univariate i-distributions. In Eq. (4), in general and e/ are the centroids of the two 
marginals of pu v -t, however a\ and 02 will approach the marginal widths <7j and a f only in the limit v — > 00, i.e. for 
the bivariate Gaussian given in Eq. (3). In-fact, the second central moments /120 = °f and /Z02 = 0/ are related to a\ 
and 02 by p2a = 7^5" CT i anc ^ = 7^2 °2 f° r v > 2. However £ remains to be the bivariate correlation coefficient. 
Exceptions to all these will occur for v < 2 and here (this happens only when A is very close to A c ) one has to use 
quartiles (i.e. spreading widths) to define ai, 02 etc.; sec [21, 22] for details. In order to test the applicability of 
the t-distribution, nuclear shell model calculations are performed for isoscalar E2 transitions in 22 Na nucleus. Fig. 1 
shows the results for A = 0.4 and 1 in the shell model hamiltonian H = h(l) + XV(2); A = 1 gives realistic nuclear 
hamiltonian. The parameters o\ and 172 in Eq. (4) are determined via their relation to p^a and /Lto2- The value of 
C = 0.88 is used as given by the exact E2 strengths. Clearly (ignoring the deviations near the ground states), the 
i-distribution gives a good description of the transition strengths with v = 6 for A = 0.4 and with a large v value, as 
expected, for A = 1. 

In larger spectroscopic spaces, instead of using a single t-distribution, to represent transition strength densities, 
it is more appropriate to partition the space. Decomposing the space into subspaces defined by h(l) eigenvalues 
£, constructing the strength distribution generated by h(l) alone, spreading this distribution by convolution with a 
t-distribution generated by V(2) and then applying some simplifying assumptions, as described in detail in [17] where 
this procedure is applied to bivariate Gaussian spreadings, it is seen that the transition strengths can be given by, 

\(E f \0\E^\ 2 = Y. M 2 {np{l-n a )) E >D{Ef) ?■ (5a) 



/+oc 
pbiv-t-,o(Ei, Ef, Si, £ f =£i-ep + e a , o\, cr 2 , C; v) dEi (5b) 
-00 

In Eq. (5a) D(Ef) denotes mean-spacing at the energy Ef, e a /3 are single particle matrix elements of O and 
(n^(l — n a )) E ' ~ {np) E% ((1 — n a )) E% with (n a ) Ei giving occupation probability for the single particle state or orbital 
a. Most remarkable is that the integral for T in Eq. (5b) can be carried out exactly for any v and this gives, 

T _ 1 

v^r(§) ^ v {o\ + 01-2^02) 



v{a\ + a\ - 2Ccrio- 2 ) 



A = E f - Ei + e - e 



(6) 
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Note that, for v > 2, o\ and 01 are related (as given above) to the marginal variances ^20 and M02 of \(Ef \ O \ Ei)\ . 
Also, the cor relation coefficient C ~ (O^VOV) / (O^O) (VV); see [15]. More importantly, as v — > 00, Eq. (6) goes 
exactly to Eq. (6) of [17] as it should be. 

To test the theory given by Eqs. (5a) and (6), numerical calculations are carried out for various A values using 25 
member EGOE(l+2) ensemble {H} = h(l) + \{V(2)} in the 924 dimensional N = 12, m = 6 space; h(l) is defined 
by the single particle energies a = + i = 1, 2, . . . , 12 and the variance of V(2) matrix elements in two-particle 

spaces is chosen to be unity. The one-body transition operator employed in the calculations is O = a\ag as in [17]. 
For the system considered, A c <~ 0.06, A^ <~ 0.2 and A t ~ 0.3. Results for six different A values, going from BW 
to Gaussian domains, are shown in Fig. 2. Clearly Eqs. (5a) and (6) obtained via the t-distribution describe the 
exact EGOE(l+2) transition strengths as we go from the BW domain with A = 0.08 to the Gaussian domain with 
A = 0.3 with v changing from 2.4 to 14; v ~ 2 — 6 for A c < A < \p and v ~ 6 — 15 for \ F < A < A t . The exact 
distributions give £ ~ 0.5 but in the fits £ is also varied (see Fig. 2) and this to some extent takes into account some 
of the approximations that led to the simple form given by Eqs. (5a) and (6). More importantly, the results in Fig. 2 
confirm that we have a good method for the calculation of transition strengths in BW domain. A calculation is also 
performed for A = 0.06 by fixing <j\ and o<i using the spreading widths of the marginals of the strength distribution 
and using ( value same as that obtained for A = 0.08. Then the deduced v value is 1.5. This and the comparisons 
in Fig. 2 clearly emphasize the role of the bivariate correlation coefficient ( in BW domain and without £ it is not 
possible to get a meaningful description (it should be mentioned that the theory in the BW domain given before 
[17, 19] uses only the marginals of the i-distribution with v = 1 and £ = 0). Thus all the problems seen before [17, 18] 
in the BW domain are cured by the bivariate-i distribution with the two parameters (f, £). 

In conclusion, random matrix ensembles generated by a mean-field plus a random two-body interaction generate 
three chaos markers. They in-turn provide a basis for statistical spectroscopy. The theory for transition strengths is 
now extended (from Gaussian domain) to BW domain down up to the A c marker by employing bivariate i-distribution. 
With atoms exhibiting a clear transition from BW to Gaussian domain (an example for Ccl to SmI atoms was shown 
in [20]), it is expected that the theory given by Eqs. (5a) and (6) will be useful in the calculation of dipole transition 
strengths in the quantum chaotic domain of atoms. 
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Shell Model ^distribution 




FIG. 1: E2 transition strengths vs {Ei, Ef). The E2 transitions considered are (0 + ,0) to (2 + ,0); with J, T and n denoting 
angular momentum, isospin and parity, nuclear levels are denoted by (J^,T). Calculations are for 6 valence nucleons in 
(2sld) shell and the hamiltonian matrix dimensions for these are 71 and 307 respectively. The hamiltonian employed is 
H{\) — h(l) + XV(2) with the single particle energies defining h(l) and two particle matrix elements defining V(2) taken from 
[23] and references therein (they define the so called Wildenthal's W-interaction) . The proton and neutron effective charges for 
the E2 operator are e p = 1.29 and e n = 0.49 respectively. All the calculations are carried out using the OXBASH computer 
code for Windows PC (2005-05 version) [24] . In the figures the energies Ei and Ef are the energies of (0 + , 0) and (2 + , 0) levels 
respectively and they are zero centered and scaled to unit width. Similarly M.E. stands for E2 transition strengths and they 
are in units of e 2 fm 4 . The vertical bars in the figures give the total strength in a given bin area; in constructing the histograms 
bin size of 0.3 is taken for both Ei and Ef. Although Eq. (4) is for EGOE(l+2), which is for spinless fermion systems, it can 
be applied directly to the shell model with good (J* ,T) states as described in [25]. 
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A = 0.08 

EGOE(l+2) v ~ 2.4, £ - 0.46 




A - 0.15 

EGOE(l+2) v = 3.6, C = 0.56 



M.E. 




A= 0.2 

EGOE(l+2) v ~ 4.5, C ~ 0.58 




FIG. 2: Transition strengths \{E f \ O \ Ei)\ 2 vs (Ei,Ef) for A = 0.08, 0.15, 0.2, 0.25, 0.28 and 0.3. In the figures Ei = E t = 
(Ei — e)/er and _B/ = _B/ = (Ef — e)/a where e and a are the centroids and widths of the state densities. Similarly M.E. stands 
for the strengths \{E f \0\ Ei)\ 2 . The EGOE(l+2) system and the one-body transition operator O are defined in the text. In 
all the calculations the strengths in the window Ei ± ^- and Ef ± ^- are summed and plotted at (Ei, Ef); A' is chosen to be 
0.1. It should be noted that the total strength is 252. As A changes from 0.08 to 0.3, the v value changes from 2.4 to 14 and 
the bivariate correlation coefficient C, changes from 0.45 to 0.62. Note the change in the scales for M.E. in the figures. 



7 



A = 0.25 

EGOE(l+2) v = 6.5, C = 0.6 



M.E. 




A = 0.28 

EGOE(l+2) v = 9, C = 0.62 



M.E. 




A =0.3 

EGOE(l+2) v = 14, C = 0.62 



M.E. 




FIG. 2 (Cont'd) 



